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ABSTRACT 

We study the rate at which gas accretes on to galaxies and haloes and investigate 
whether the accreted gas was shocked to high temperatures before reaching a galaxy. 
For this purpose we use a suite of large cosmological, hydro dynamical simulations from 
the OWLS project, which uses a modified version of the smoothed particle hydrody- 
namics code gadget- 3. We improve on previous work by considering a wider range 
of halo masses and redshifts, by distinguishing accretion on to haloes and galaxies, by 
including important feedback processes, and by comparing simulations with different 
physics. 

Gas accretion is mostly smooth, with mergers only becoming important for groups 
and clusters. The specific rate of gas accretion on to haloes is, like that for dark 
matter, only weakly dependent on halo mass. For halo masses Mh a io ^> 10 11 M it is 
relatively insensitive to feedback processes. In contrast, accretion rates on to galaxies 
are determined by radiative cooling and by outflows driven by supernovae and active 
galactic nuclei. Galactic winds increase the halo mass at which the central galaxies 
grow the fastest by about two orders of magnitude to Mh a io ~ 10 12 M . 

Gas accretion is bimodal, with maximum past temperatures either of order the 
virial temperature or < 10 5 K. The fraction of gas accreted on to haloes in the hot 
mode is insensitive to feedback and metal-line cooling. It increases with decreasing 
redshift, but is mostly determined by halo mass, increasing gradually from less than 
10% for rsj 10 11 M to greater than 90% at ~ 10 13 M . In contrast, for accretion on to 
galaxies the cold mode is always significant and the relative contributions of the two 
accretion modes are more sensitive to feedback and metal-line cooling. On average, 
the majority of stars present in any mass halo at any redshift were formed from gas 
accreted in the cold mode, although the hot mode contributes typically over 10% for 
M halo > 10 11 M . 

Thus, while gas accretion on to haloes can be robustly predicted, the rate of 
accretion on to galaxies is sensitive to uncertain feedback processes. Nevertheless, it 
is clear that galaxies, but not necessarily their gaseous haloes, are predominantly fed 
by gas that did not experience an accretion shock when it entered the host halo. 

Key words: galaxies: evolution - galaxies: formation - intergalactic medium - cos- 
mology: theory 



1 INTRODUCTION 

In the standard cosmological constant or vacuum dominated 
cold dark matter (ACDM) model mass assembles hierarchi- 
cally, with the smallest structures forming first. While the 
collapse of dark matter halts as it reaches virial equilibrium 
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in haloes, baryons can radiate away their binding energy, 
allowing them to collapse further and fragment into smaller 
structures, such as stars and galaxies. These galaxies then 
grow through mergers and gas accretion. 

After virialization, a gas cloud can be in one of three 
regimes. If the characteristic cooling time exceeds the Hub- 
ble time-scale, the gas will not be able to radiate away the 
thermal energy that supports it and will therefore not col- 



© 2010 RAS 



2 F. van de Voort et al. 



lapse. If, on the other hand, the cooling time is smaller 
than the Hubble time, but larger than the dynamical time- 
scale, then the cloud can adjust its density and tempera- 
ture quasi- statically. It will increase both its density and 
temperature while maintaining hydrostatic equilibrium. Fi- 
nally, if the cooling time is shorter than the dynamical time, 
the cloud will cool faster than it can collapse, lowering the 
Jeans mass and possibly leading to fragmentatio n. This is 
the regime in which galaxies are thought to form (|Silldll977l : 
In reality the situation must, how- 
ever, be more complicated as the density, and thus the cool- 
ing and dynamical times, will vary with radius. 

Gas falling towards a galaxy gains kinetic energy until 
it reaches the hydrostatic halo. If the infall velocity is su- 
personic, it will experience a shock and heat to the virial 
temperature of the halo. According to the simplest pic- 
ture of spherical collapse, all gas in a dark matter halo is 
heated to the virial temperature of that halo, reaching a 
quasi-static equilibrium supported by the pressure of the 
hot gas. Gas can subsequently cool radiatively and settle 
into a rotationally supported disc, where it can form stars 
(e.g. iFall fe Efstathioul ll980h. We call this form of gas ac- 
cretion 'hot accretion' (|Katz et al.ll2003l ; IKeres et al.ll2005h . 

Within some radius, the so-called cooling radius, the 
cooling time of the gas will, however, be shorter than the 
age of the Universe. If the cooling radius lies well inside the 
halo, which is the case for high-mass haloes, a quasi-static, 
hot atmosphere will form. Accretion on to the galaxy is then 
regulated by the cooling function. If, on the other hand, this 
radius is larger than the virial radius, then there will be no 
hot halo and the gas will not go through an accretion shock 
at the virial radius. Because gas accreted in this manner 
may never have been heated to the virial temper ature, we 
refer t o this mode of acc retion as 'cold accretion' (|Katz et al.l 
120031 : IKeres et al]l2005T ). The rapi d cooling of gas i n low - 
mass haloes was already shown by Rees & Ostrikerl (|l977T ) 
and I White fc Reed (|l978T ). The accretion rate on to the cen- 
tral galaxy t hen depends on the in fall rate, but not on the 
cooling rate (|White fe Fr enk 1991). Simulations confirmed 
the existence of gas inside haloes which was never heated 
to the virial t emperatu re of the halo (iKatz &; Gun n 1991; 
iKav et alJlipOol ; iFardal et al1l200ll ; IKatz et al.ll2003l ). 

The cooling time t coo i is shorter at higher redshift, 
because the density, p, is higher and t coo \ oc The 
Hubble time £h is also shorter, but with a weaker depen- 
dence, £h oc p -1 ^ 2 . Hence, the mode of gas accretion on 
to galaxies will depend on redshift, with cold accretion 
more prevalent at higher redshifts for a fixed virial tem- 
perature. The mode of gas accretion will also depend on 
halo mass. Higher-mass haloes have higher virial tempera- 
tures and thus, at least for T > 10 6 K, longer cooling times. 
Hence, hot accretion will dominate for the most massive 
haloes. The dominant form of accretion will thus depend 
on b oth the mass of the halo and on the accretion red- 
shift (jKatz et al.ll2003l;lBirnboim fc Dekel 20o |;lKeres et al 
20051 : lOcvirk et al.1 120081 : IKeres et all l2009al : iBrooks et al 
2009l : ICrain et alJl20ld ). 



Although both analytic and semi-analytic studies 
of g alaxy format i on usually assume spherical symmetry 
(e.g lBinnevlll977l : IWhite fe Frenklll99ll : iBirnboim & Dekell 
2003), numerical simulations show different geometries. As 
the matter in the Universe collapses, it forms a network 



of sheets and filaments, the so-called 'cosmic web'. Galax- 
ies form in the densest regions, the most massive galaxies 
where filaments intersect. These structures can have an im- 
portant effect on gas accretion. If a galaxy is being fed along 
filaments, the average density of the accreting gas will be 
higher. The cooling time will thus be smaller and it will be 
easier to radiate away the gravitational binding energy. Fil- 
aments may therefore feed galaxies preferentially through 
cold accretion. Both modes can coexist. Especially at high 
redshift, cold str eams penetrate the hot virialized haloes of 
massive galaxies (iKeres et a l. 2005; Dekel & Birnboim 2006; 
lOcvirk et aT1l2008l : IKeres et al.ll2009al : iDekel et al.ll2009ah . 

Galaxy colours and morphologies are observed to be 
roughly bimodal. They can be divided into two main classes: 
blue, star forming spirals and red, passiv e ellipticals (e.g. 
iKauffmann et al.ll2003l : iBaldrv et al.ll2004h . The latter tend 
to be more massive and to reside in c lusters. It has been 
suggested by IDekel &; Birnboiml (|2006l ) that this observed 
bimodality is caused by the two different mechanisms for 
gas accretion. If, for example, feedback from active galac- 
tic nuclei (AGN) were to prevent the hot gas from cooling, 
massive galaxies with little cold accretion would have very 
low star formation rates. If, on the other hand, cold accre- 
tion flows were less susceptible to such feedback, then low- 
mass haloes could host discs of cold gas, which could form 
stars efficiently. However, this cannot be the whole story, 
because suppressing star formation from gas accreted in the 
hot m ode is not nearly s ufficient to reproduce the observa- 
tions (jKeres et al.ll2009bh . 

In this work we use a large number of cosmological, hy- 
drodynamical simulations fro m the Overwhelm ingly Large 
Simulations project fOWLS: ISchave et al.ll2010l ) to investi- 
gate accretion rates, the separation of hot and cold modes, 
and their dependence on halo mass and redshift. Our work 
extends earlier work in several ways. By combining simula- 
tions with different box sizes, each of which uses at least as 
many particles as previous simulations, we are able to cover a 
large dynamic range in halo masses and a large range of red- 
shifts. While earlier work used only a single physical model, 
our use of a range of models allows us to study the role of 
metal-line cooling, supernova (SN) feedback, and AGN feed- 
back. Some of these p rocesses had been i gnored by earlier 
studies. For example, IKeres et all (|2005l . l2009al lry) ignored 
metal-l ine cooling and SN feed back. Using a semi-analytic 
model, iBenson &; Bowerl (|2010l ) em phasized the import ance 
of including the effect of feedback. IBrooks et al.l (|2009h did 
not include metal- line cooling, but did include SN feedback. 
lOcvirk et al.l (|2008l ) included both metal-line cooling and 
weak SN feedback, but could only study the high-redshift 
behaviour as their simulation was stopped at z ~ 1.5. AGN 
feedback was not inclu ded by any previous study, although 
iKhalatvan et al.l |2008) did show the existence of cold ac- 
cretion in a simulation of a single group of galaxies with 
AGN feedback. In contrast to previous work, we distinguish 
between accretion on to haloes and accretion on to galax- 
ies (i.e. the interstellar medium). We find that while hot 
mode accretion dominates the growth of high-mass haloes, 
cold mode accretion is still most important for the growth 
of the galaxies in these haloes. The different physical mod- 
els give similar results for accretion on to haloes, implying 
that the results are insensitive to feedback processes, but 
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Table 1. Simulation parameters: simulation identifier, comoving box size (Lbox)? number of dark matter particles (N, the number of 
baryonic particles is equal to the number of dark matter particles), mass of dark matter particles (?tidm)) initial mass of gas particles 
(m g as), final simulation redshift, number of resolved haloes (containing more than 250 dark matter particles) at z = 2 (Ni ia \ (z = 2)), 
and number of resolved haloes (containing more than 250 dark matter particles) at z = (iVh a io(2 = 0)). 



simulation 




N 


WDM 




Zfinal 


^halo(£ = 2) 


A^halo^ = 0) 




(h^Mpc) 




(M ) 


(M ) 








L100N512 


100 


512 3 


5.56 x 10 8 


1.19 x 10 8 





10913 


18463 


L WON 256 


100 


256 3 


4.45 x 10 9 


9.84 x 10 9 





838 


2849 


L050N512 


50 


512 3 


6.95 x 10 7 


1.48 x 10 7 





12648 


15884 


L025N512 


25 


512 3 


8.68 x 10 6 


1.85 x 10 6 


2 


12768 




L025N256 


25 


256 3 


6.95 x 10 7 


1.48 x 10 7 


2 


1653 




L025N128 


25 


128 3 


5.56 x 10 8 


1.19 x 10 8 





159 


338 



the inclusion of such processes is important for accretion on 
to galaxies. 

This paper is organized as follows. The simulations are 
described in Section [2] including all the model variations. 
In Section [3] we describe some properties of the gas in the 
temperature-density plane for our fiducial simulation. The 
ways in which haloes are identified and gas accretion is de- 
termined is discussed in Section [4] The accretion rates for 
all simulations can be found in Section [5j for accretion on 
to haloes as well as for accretion on to galaxies. In Section [6] 
we compare hot and cold accretion on to haloes for the dif- 
ferent simulations and we do the same for accretion on to 
galaxies in Section Finally, we compare our results with 
previous work in Section [8] and summarize our conclusions 
in Section [9] 

2 SIMULATIONS 

To investigate the temperature distribution of accreted 
gas, we u se a modified version of GADGET-3 (last de- 
scribed in Springel 2005), a smoothed particle hydrody- 
nami cs (SPH) code that uses t he entropy formulation of 
SPH ([Springel &; Hernquistll 20021 ), which conserves both en- 
ergy and e ntropy where appro priate. This work is part of 
the OWLS (|Schave et al.ll2010h project, which consists of a 
large number of cosmological simulations, with varying (sub- 
grid) physics. We first summarize the subgrid prescriptions 
for the reference simulation. The other simulations are de- 
sc ribed in Section 12. 1[ As the simulations are fully described 
in ISchave et al.l ([2010), we will only summarize their main 
properties here. 

The cosmological simulations described here assume a 
ACDM cosmology with parameters as determined from the 
Wilkinson microwave anisotropy probe year 3 (WMAP3) 
results, On = 1 - Ov = 0.238, Q h 0.0418, h 0.73, 
as = 0.74, n — 0.951. These values are con sistent with the 
WMAP year 7 data (|Komatsu et alJl201lh . 

A cubic volume with periodic boundary conditions 
is defined, within which the mass is distributed over 
N 3 dark matter and as many gas particles. The box 
size (i.e. the length of a side of the simulation vol- 
ume) of the simulations used in this work is 25, 50, 

1 The only significant discrepancy is in as, which is 8% lower 
than the value favoured by the WMAP 7-year data. 



or 100 h x Mpc, with N — 512, unless stated other- 
wise. The (initial) particle masses for baryons and dark 
matter are 1.2 x 10 8 ( — L h b -T Mpc ) 3 ( ^)" 3 M and 5.6 x 

1q8 ( ioo t-°i X Mpc ) 3 (5ll)~ 3 M ©' respectively, and are listed in 
Table [1] as are the number of resolved haloes at z — 2 
and z = 0. We use the notation L***N###, where *** 
indicates the box size and ### the number of parti- 
cles per dimension. The gravitational softening length is 

7 - 8 ( iooh-?M P c )(^)" 1 ft " lk P c comoving, i.e. 1/25 of the 
mean dark matter particle separation, but we imposed a 
maximum of 2 ( — t-T Mpc )( sf^)" 1 /i _1 kpc proper. 

The primordial abundances are X = 0.752 and Y = 
0.248, where X and Y are the mass fractions of hydrogen and 
helium, respectively. The abundances of eleven elements (hy- 
drogen, helium, carbon, nitrogen, oxygen, neon, magnesium, 
silicon, sulphur, calcium, and iron) released by massive stars 
(type II SNe and stellar winds) and intermediate mass stars 
(type la SNe a nd asymptotic giant bra nch stars) are followed 
as described in lWiersma et al.l (12009 q). We a ssum e the stel- 
lar initial mass function (IMF) of Chabrier (2003), ranging 
from 0.1 to 100 M . As described in lWiersma et alJ (|2009ah , 
radiative cooling and heating are computed element-by- 
element in the pres ence of the cosmic microw ave background 
radiation and the lHaardt &; Madaul (|200ll ) model for the 
UV/X-ray background from galaxies and quasars. 

Star formation is modelle d according to the recipe of 
ISchave fe Dalla Vecchial (|2008h . The Jeans mass cannot be 
resolved in the cold, interstellar med ium (ISM), which coul d 
lead to artificial fragmentation (e.g. Bat e &; Burkertlfl997l ). 
Therefore a polytropic equation of state P to t oc Pga? is im- 
plemented for densities exceeding nu = 0.1 cm -3 , where 
Ptot is the total pressure and p ga s the density of the gas. It 
keeps the Jeans mass fixed with respect to the gas density, as 
well as the ratio of the Jeans length and the SPH smoothing 
kernel. Gas particles with proper densities nu ^ 0.1 cm -3 
and temperatures T ^ 10 5 K are moved on to this equation 
of state and can be converted into star particles. The star 
formation rate (SFR) per unit mass depends on the gas pres- 
sure and is set to repr oduce the observed Kennicutt-Schmidt 
law (|Kennicuttlll998l ). 

Feedback fr om star formation is imp l ement ed using the 
prescription of iDalla Vecchia k, Schavd (|2008l ). About 40 
percent of the energy released by type II SNe is injected lo- 
cally in kinetic form. The rest of the energy is assumed to be 
lost radiatively. Each gas particle within the SPH smoothing 
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Table 2. Simulation parameters: simulation identifier, cooling including metals (Z cool), wind velocity (f w i n d); wind mass loading (77), 
and AGN feedback included (AGN). Differences from the reference model are indicated in bold face. The last column gives the box size 
and particle number, from Tabled used with these simulation parameters (Lbox an< ^ 



simulation 


Z cool 


^wind 


77 


AGN 


L box and N 






(km/s) 








REF 


yes 


600 


2 


no 


all listed in Tabled 


NOSN 


yes 








no 


L100N512 at z = 2 


NOSN.NOZCOOL 


no 








no 


L100N512, L025N512 


NO Z COOL 


no 


600 


2 


no 


L100N512, L025N512 


WML4 


yes 


600 


4 


no 


L100N512, L025N512 


WML1V848 


yes 


848 


1 


no 


L100N512, L025N512 


DBLIMF- 


yes 


600 & 1618 


2 


no 


L100N512, L025N512 


CONTSFV1618 












WDENS 


yes 


density dependent 


no 


L WON 5 12, L025N512 


AGN 


yes 


600 


2 


yes 


L WON 512, L025N512 



kernel of the newly formed star particle has a probability of 
being kicked. For the reference model, the mass loading pa- 
rameter 77 = 2, meaning that on average the total mass of the 
particles being kicked is twice the mass of the star particle 
formed. Because the winds sweep up surrounding material, 
the effective mass loading can be higher. The initial wind 
velocity is 600 km/s for the reference model, which is con- 
sistent with observat ions of starburst galaxie s, both locally 
(IVeilleux et al.1 12005) and at high redshift (Sha plev et al.l 
2003). ISchave et al.l (|2010h showed that these parameter val- 
ues yield a peak, global star formation rate density that 
agrees with observations. 

The simulation data is saved at discrete output redshifts 
with interval Az = 0.125 at < z ^ 0.5, Az = 0.25 at 
0.5 < z < 4, and Az = 0.5 at 4 < z < 9. This is the time 
resolution used for determining accretion rates. 



2.1 Model variations 

To see whether or not our results are sensitive to specific 
physical processes or subgrid prescriptions, we have per- 
formed a suite of simulations in which many of the simu- 
lation parameters are varied. These are listed in Table [2] 

The importance of metal-line cooling can be demon- 
strated by comparing the reference simulation (REF) to 
a simulation in which primordial abundances are assumed 
when calculating the cooling rates (NOZCOOL). Similarly, 
the effect of including SN feedback can be studied by com- 
paring a simulation without SN feedback (NOSN) to the ref- 
erence model. Because the metals cannot be expelled with- 
out feedback, they pile up, causing efficient cooling and 
star formation. To limit the cooling rates, we performed 
a simulation in which both cooling by metals and feed- 
back from SNe were ignored (NOSN.NOZCOOL). Simula- 
tion NOSN.L025N512 was stopped just below z = 4 and 
the z = 0.125 snapshot is missing for NOSN.L100N512. We 
therefore cannot always show results for this simulation. 

In massive haloes the pressure of the ISM is too high 
for winds with velocities of 600 km/s to blo w the gas out 
of the galaxy (|Dalla Vecchia &; Schavd 12008). Keeping the 
wind energy per unit stellar mass constant, we increased 
the wind velocity by y/2 to v w ind = 848 km/s, while halving 
the mass loading to 77 = 1 (WML1V848). To enable the 



winds to eject gas from haloes with even higher masses, the 
velocity and mass loading can be scaled with the local sound 
speed, while keeping the energy injected per unit stellar mass 
constant (WDENS). 

To investigate the dependence on SN energy given to 
the ISM, we ran a simulation using almost all of the energy 
available from SNe. The wind mass loading 77 = 4, a factor 
of two higher than in the reference simulation. The wind 
velocity is the same, v w ind = 600 km/s (WML4). 

There is some evidence that the IMF is top- 
heavy in e xtreme enviro n ments , like starburst galax- 
ies (e.g. IPadoan et al.l Il997l ; iKlessen et al.1 120071 : 
iDabringhausen et al.l [2009). We have performed a sim- 
ulation in which the IMF is top-heavy for stars formed at 
high pressures, Ptot/ku > 1.7 x 10 5 cm -3 K, where ku is 
Boltzmann's constant, and therefore in high-mass haloes. 
The critical pressure was chosen such that ~10% of the 
stars are formed with a top-heavy IMF. It results in more 
SNe per solar mass of stars formed and therefore stronger 
winds. More energy is therefore put into the winds blown 
by these star particles, v w ind = 1618 km/s and 77 = 2 
( DBLIMFCONTSF VI 6 18). 

Finally, we have included AGN feedback (AGN). Black 
holes inject 1.5% of the rest-mass energy of the accreted gas 
into the surrounding matt er in the form of heat. The model 
is described and tested in iBooth &; Schavd (|2009r ) , who also 
demonstrate that it reproduces the observed mass density in 
black holes and the observed scaling relations between black 
hole mass and bo th central stellar veloc ity dispersion and 
and stellar mass. iMcCarthv et al . (2010) have shown that 
our model AGN reproduces the observed stellar mass frac- 
tions, SFRs, stellar age distributions and thermodynamic 
properties of galaxy groups. 



2.2 Maximum past temperature 

The Lagrangian nature of the simulation is exploited by 
tracing each fluid element back in time, which is particu- 
larly convenient for this project, in which we are studying 
the temperature history of accreted gas. During the simu- 
lations the maximum temperature, T max , and the redshift 
at which it was reached, z max , were stored in separate vari- 
ables. The variables were updated for each SPH particle at 
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Figure 1. Current temperature (left) and maximum past temperature (right) against current overdensity for gas particles at z = 2 
in the reference simulation in a 50 h~ 1 Mpc box. The logarithm of the total gas mass in a pixel is used for colour coding. The black 
lines separate the different phases of the gas. Gas with T > 10 5 K has been shock heated. The T — p relation for cold, low-density gas 
(T < 10 5 K and p < 10(p) is set by heating by the UV background and adiabatic cooling. At p/{p) > 10 1-5 radiative cooling dominates 
over cooling by the expansion of the Universe. Gas with tih > 0.1 cm -3 (p/ (p) > 10 4 3 at z = 2) is assumed to be part of the unresolved, 
multiphase ISM and is put on an effective equation of state. The 'temperature' of this gas merely reflects the pressure of the multiphase 
ISM and is therefore not used to update T ma x- The scatter in the equation of state is caused by adiabatic cooling of inactive particles in 
between time steps. Gas that is condensing on to haloes or inside haloes has T max > 10 4 5 K, which reflects the peak in the T-p relation 
(at p/(p) > 10 15 , T ~ 10 4 - 5 K) visible in the left panel. 



every time step for which the temperature was higher than 
the previous maximum temperature. The artificial temper- 
ature the particles obtain when they are on the equation 
of state (i.e. when they are part of the unresolved multi- 
phase ISM) was ignored in this process. This may cause us 
to underestimate the maximum past temperature of gas that 
experienced an accretion shock at higher densities. Ignoring 
such shocks is, however, consistent with our aims, as we are 
interested in the maximum temperature reached before the 
gas entered the galaxy. 

Shock heating could be missed if the particle heats 
and cools rapidly in a single time step. In SPH simula- 
tions, a shock is smeared out over a few smoothing lengths , 
leading to in-shock cooling (jHutchings fe Thomas! (2000). 
We find the average value for the infall velocity to be 
~ 5 x 102 ( lo i3 M M0 ) 1/3 (l + z) 1/2 km/s. On average the 
SPH smoothing length at halo accretion is ~ 10 2 kpc (co- 
moving) for all L100N512 simulations and at all redshifts. 
The particle will take 10 8 ( 10 i a M M )~ 1/3 (l + a)~ 3/2 years to 
traverse this distance. A shock will therefore take at least a 
few times 10 8 years. In reality the accretion shock will pro- 
ceed almost instantaneously, minimizing radiative losses. In 
the simulations described here gas cools at each time step, 
so also while it is being shocked. The finite spatial resolution 
can thus result in lower maximum temperatures than the ac- 
tual post-shock temperatures. If this effect were important, 
increasing the resolution would increase the hot fraction, 
i.e. the fraction of the gas that is accreted in the hot mode, 
because the accretion shocks would be less broadened. The 
opposite is the case, as we will show in Figure [9] 

Even with infinite time resolution, the post-shock tem- 
peratures may, however, not be well defined. Electrons and 
protons will temporarily have different temperatures in the 
post-shock gas, because they differ in mass, an effect that 
we have not included. It will take a short while before they 



equilibrate through collisions or plasma effects. Another ef- 
fect, which is also not included in our simulations, is that 
shocks may be preceded by radiation from the shock, which 
may affect the temperature evolution. 



3 THE TEMPERATURE-DENSITY 
DISTRIBUTION 

The left-hand panel of Figure [1] shows the mass- weighted dis- 
tribution of gas in the temperature-density plane at z — 2 for 
simulation REF_L050N512. Gas with densities up to ~ 10 2 
times the cosmic average density represents the diffuse in- 
tergalactic medium (IGM) . A significant fraction of this gas 
resides in filamentary structures. It can be heated to tem- 
peratures above 10 5 K when kinetic energy, generated by 
gravitational infall or galactic winds, is converted into ther- 
mal energy. We refer to this tenuous, shock-heated gas as 
the warm- hot intergalactic medium (WHIM). The intraclus- 
ter medium (ICM) is the very hot T > 10 7 K gas located in 
galaxy groups and clusters. Gas at overdensities p/ (p) > 10 2 , 
but much lower temperatures (T ~ 10 4 K) resides mostly in 
filaments and low- mass haloes. 

Most of the gas is located in the purple region, with 
p/{p) < 10. The temperature of this gas is determined by 
the combination of photoheating by the UV background and 
adiabatic cooling by the expansion of the Universe. Although 
the slope of this temperature-density relation is close to adi- 
abatic, it is actually determined by the temperature d epen- 
dence of the recombination rate (|Hui fe Gnedin|[l997h . 

The turnover density, above which the typical gas tem- 
peratures decrease, occurs when radiative cooling starts to 
dominate over adiabatic cooling. The overdensity at which 
this happens depends on redshift. At z = 2 it occurs at 
p/{p) « 10 15 . The distribution of the WHIM (broad red 
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region, with T > 10 5 K) is set in part by the cooling func- 
tion, to which especially heavier atoms, like oxygen, con- 
tribute. In particular, the lack of dense gas (p/{p) > 10 3 ) 
with 10 5 K < T < 10 7 K is due to radiative cooling. 

Gas with proper hydrogen number density 71h > 
0.1 cm -3 , corresponding at this redshift to overdensity 
P/{p) > 10 4 ' 3 , represents the ISM. This high-density gas 
is put on an equation of state if its temperature was below 
10 5 K when it crossed the density threshold, because the cold 
and warm phases of this dense medium are not resolved by 
the simulations. Therefore, the temperature merely reflects 
the imposed effective pressure and the density should be in- 
terpreted as the mean density of the unresolved multi-phase 
ISM. The spread in the temperature-density relation on the 
equation of state is caused by the adiabatic extrapolation 
of inactive particles between time steps. In addition, the 
relation is broadened by differences in the mean molecular 
weight, /i, of the gas, which depends on the density, temper- 
ature and elemental abundances of the gas. 

The right-hand panel of Figure [1] shows the maximum 
past temperature reached at z ^ 2, as a function of the 
z — 2 baryonic overdensity. All dense gas has reached tem- 
peratures of > 10 4-5 K at some point. Dense (p/{p) > 10 1 " 5 ) 
gas cannot have T max much below 10 4-5 K, because of pho- 
toheatin£|. The maximum past temperature tends to re- 
main constant once the gas has reached densities > 10 2 (p), 
resulting in the horizontal trend in the figure. There is no 
dense p/{p) ~ 10 4 gas at the highest maximum temperatures 
(above 10 7 K), because the cooling time of gas that is heated 
to this temperature at lower densities (p/(p) ~ 10 2-3 ) is 
longer than a Hubble time, preventing it from cooling and 
reaching higher densities. 

The maximum past temperature reached for dense gas 
depends on whether the gas has been heated to the virial 
temperature of its halo and on whether it has been shock- 
heated by galactic winds. 



4 DEFINING GAS ACCRETION 

To see how haloes accrete gas, we first need to find and select 
the haloes. This can be done in several different ways and we 
will discuss three of them. Although we choose to use the 
one based on the gravitational binding energy, our results 
are insensitive to the halo definition we use. Finally, we link 
haloes in two subsequent snapshots in order to determine 
which gas has entered the haloes. 

4.1 Identifying haloes and galaxies 

The first step towards finding gravitationally bound struc- 
tures is to identify dark matter haloes. These can for exam- 
ple be found using a Friends-of-Friends (FoF) algorithm. If 
the separation between two dark matter particles is less than 
20% of the average separation (the linking length b = 0.2), 
they are placed in the same group. Because the particles all 
have the same mass, a fixed value of b will correspond to a 

2 Exceptions are gas that reached high densities before reioniza- 
tion, which happens at z = 9 in our simulations, and gas with 
very high metallicities. 



fixed overdensity at the boundary of the group of p/ (p) w 60 
(e.g. Frenk et al. 1988). Assuming a radial density profile 
p(r) oc r -2 , corresponding to a flat rotation curve, such a 
group has an avera ge overdensity of (phaio) / (p) ~ 180 (e.g. 
lLacev fc C ole 1994), close to the value for a virialized ob- 
ject predicted by th e top-hat spherical collapse model (e.g. 
[Padmanabhanl 2002). The minimum number of dark matter 
particles in a FoF group is set to 25. Many of the smallest 
groups will not be gravitationally bound. Baryonic particles 
are placed in a group if their nearest dark matter neighbour 
is part of the group. 

Problems arise because unbound particles can be at- 
tached to a group, physically distinct groups can be linked 
by a small (random) particle bridge, and because substruc- 
ture within a FoF halo is not identifie d. We use subfind 
(jSpringel et al.ll200ll : iDolag et al1l2009h on the FoF output 
to find the gravitationally bound particles and to identify 
subhaloes. The properties of gas that is being accreted are 
expected to depend on the properties of the parent (or main) 
halo in which the subhaloes are embedded. We therefore only 
look at accretion on to these main haloes and exclude gas 
accretion by satellites. 

Another way of defining haloes is to use a spherical 
overdensity criterion. The radius R vir , centred on the most 
bound particle of a FoF halo, is found within which the aver- 
age density agrees with the prediction of t he top-hat spheri- 
cal co llapse model in a ACDM cosmology ([Bryan Sz Normanl 
1998). All the particles within R vir are then included in the 
halo. 

We chose to use only main haloes identified by sub- 
find, but we have checked that our results do not change 
significantly when using FoF groups or spherical overdensi- 
ties instead. 

Except for Figures [2] [9] and 1151 we include only main 
haloes containing more than 250 dark matter particles in our 
analysis of halo accretion. This corresponds to a minimum 
total halo mass of M ha i « 10 112 M in the 100 /i _1 Mpc 
box, 10 10 - 3 M in the 50 /i _1 Mpc box, and 10 9 4 M in the 
25 /i -1 Mpc box. For these limits our mas s functions agree 
very well with the lSheth fe Tormenl (|l999h fit. 

For each resolved halo, we identify the ISM of the cen- 
tral galaxy with the star forming (i.e. nu > 0.1 cm -3 ) gas 
particl es in the main hal o which are inside 15% of the virial 
radius (Sales et al. 2010|). We use 15% of the virial radius to 
exclude small star forming substructures in the outer halo, 
which are not identified by subfind. Gas accretion on to 
satellite galaxies is excluded. 

Because the galaxy is much smaller than its parent halo, 
it is not as well resolved. When investigating accretion on 
to galaxies (as opposed to haloes), we therefore impose a 
1000 dark matter particle limit, corresponding to a minimum 
total halo mass of M ha i « 10 118 M in the 100 /i _1 Mpc 
box, 10 10 - 9 M in the 50 fc _1 Mpc box, and 10 10 M in the 
25 /i _1 Mpc box. The exceptions are Figure [2l and ITol which 
show accretion rates and hot fractions down to 1 dex below 
this limit. See Section [5] for more discussion on convergence. 

4.2 Selecting gas particles accreted on to haloes 

We select gas particles accreted on to haloes as follows. For 
each halo at z = Z2 (which we will also refer to as 'the de- 
scendant') we identify its progenitor at the previous output 
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redshift z\ > 22. We determine which halo contains most 
of the descendant's 25 most bound dark matter particles 
and refer to this halo as 'the progenitor'. If the fraction of 
the descendant's 25 most bound particles that was not in 
any halo at z\ is greater than the fraction that was part of 
the progenitor, then we discard the halo from our analysis, 
which rarely happens above our resolution threshold. If two 
or more haloes contain the same number of those 25 parti- 
cles, we select the one that contains the dark matter particle 
that is most bound to the descendant. 

We identify those particles that are in the descendant, 
but not in its progenitor as having been accreted at Z2 ^ z < 
z\. The accreted particles have to be gaseous at zi, i.e. before 
they were accreted, but can be either gaseous or stellar at 
redshift 22. The accreted gas can have densities exceeding 
the star formation threshold, in which case it cannot obtain 
a higher maximum past temperature. Gas can be accreted 
multiple times. 

To distinguish mergers from smooth accretion, we ex- 
clude accreted particles that reside in (sub) haloes above 
some maximum mass at redshift z\ . We would like to have a 
criterion that is not directly dependent on resolution, so that 
the same objects are included in runs with different particle 
numbers if the simulations are converged with respect to res- 
olution. We therefore set this maximum allowed halo mass 
to 10% of the descendant's (main halo) mass. Thus, smooth 
accretion excludes mergers with a mass ratio greater than 
1:10. We experimented with different thresholds and the re- 
sults are insensitive to this choice. 



4.3 Selecting gas particles accreted on to galaxies 

We consider particles that are part of the ISM or stellar at 
22, and that were gaseous but not part of the ISM at zi, 
to have been accreted on to a galaxy at Z2 ^ z < z\. Gas 
can be accreted multiple times. Accretion on to the ISM and 
accretion on to a galaxy are the same in this study and these 
terms are used interchangeably. 

In this way, accretion through galaxy mergers is auto- 
matically excluded, because the gas that is part of the ISM 
of another galaxy at z\ is excluded. We can allow for merg- 
ers by identifying those particles that are in the galaxy at 
Z2 , but were not in the progenitor galaxy at z\ . 



5 TOTAL GAS ACCRETION RATES 

For brevity, accretion between, for example, z — 2.25 and 
z — 2 will be referred to as accretion just before z — 2. 

The average total gas accretion rate of a halo depends 
on its mass. Figure [2] shows the average specific accretion 
rates QmMgas/^bMhaio, where Q m and Qb are the matter 
and baryon density parameters, respectively, as a function 
of halo mass for z — 2 (top panel) and z — (bottom 
panel) for haloes containing at least 100 dark matter parti- 
cles for various simulations using the reference model. Thick 
curves show the specific accretion rates on to haloes. We 
have divided the gas accretion rate by the total halo mass 
and baryon fraction / , so that the normalization gives 
an estimate of the time it would take a halo to grow to its 
current mass with the current accretion rate. 
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Figure 2. Specific gas smooth accretion rates on to haloes 
(higher, thick curves) and central galaxies (lower, thin curves) 
against total halo mass at z = 2 (top panel) and z = (bottom 
panel). The normalization for halo accretion gives an estimate 
of the time it would take a halo to grow to its current mass at 
the current accretion rate. The same is not true for galaxy accre- 
tion, since we use the same normalization as for halo accretion. 
Both panels show simulations spanning a factor 64 in mass res- 
olution. Each mass bin contains at least 10 haloes. Big (small) 
arrows correspond to the adopted resolution limit for accretion 
on to haloes (galaxies) for L025N512 (red), L050N512 (blue), 
L WON 512 (black) and L WON 256 (green). The specific halo ac- 
cretion rate is converged and increases slightly with halo mass. 
The specific galaxy accretion rate is not fully converged at z = 2, 
but the convergence is better at z = 0. The galaxy accretion rate 
increases with halo mass for Mh a i < 10 12 Mq and decreases 
with halo mass for higher halo masses. It is much smaller than 
the halo accretion rate. 



The thin curves in Figure [2] show specific smooth ac- 
cretion rates on to galaxies against halo mass. While the 
inverse of the specific halo accretion rate equals the time it 
takes to grow the halo at its current accretion rate, the same 
is not true for the specific galaxy accretion rate, because we 
divide by the halo mass rather than by the galaxy mass. 
This definition allows us, however, to directly compare the 
two accretion rates. 

Comparing the three thick curves in each panel, we see 
that the specific halo accretion rate is converged with the 
numerical resolution. It increases with halo mass, especially 
at low halo masses. The specific galaxy accretion rate (thin 
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curves) is not fully converged at z = 2, but the conver- 
gence improves at z = (compare the two highest resolu- 
tions, dotted and solid curve, at the high mass end). Below 
Mhaio < 10 12 M© the specific galaxy accretion rate increases 
more steeply with halo mass than the specific halo accre- 
tion rate. Above this halo mass, the specific galaxy accre- 
tion rate decreases steeply, whereas the rate keeps increasing 
for haloes. The accretion rate on to galaxies is always much 
lower than the halo accretion rate and this difference is larger 
at z = 0. Hence, only a small fraction of the gas accreted 
on to haloes ends up in galaxies and the efficiency of galaxy 
formation is highest in haloes with Mhaio ~ 10 12 M©. 

The dynamic range covered by these simulations is very 
large, because we use different box sizes. The resolution in- 
creases with decreasing box size. We have performed box 
size convergence tests (at fixed resolution), but do not show 
them as the convergence with box size is excellent for all 
halo masses. Increasing the resolution shows that the halo 
accretion rates are not fully converged around the 100 parti- 
cle resolution threshold adopted in Figure [2] In this regime, 
halo mergers cannot always be identified. In the rest of this 
paper, we therefore set the minimum halo mass for accretion 
on to haloes to correspond to 250 dark matter particles, as 
indicated by the big arrows. 

The galaxy accretion rates diverge at the low-mass end. 
This is expected, because haloes with 100 dark matter parti- 
cles will have very few star forming gas particles and because 
15% of the virial radius is close to the spatial resolution 
limit. We will therefore only include the haloes where 15% 
of the virial radius is larger than 5 times the softening length 
and therefore consisting of at least 1000 dark matter parti- 
cles, as indicated by the small arrows. The galaxy accretion 
rates at the high-mass end are not completely converged at 
z = 2. The accretion rate decreases by up to a factor of two 
if the mass resolution is increased by a factor of 64. With in- 
creasing resolution, more gas reaches higher densities. If this 
gas becomes star forming before accreting on to the central 
galaxy, it is not included in the smooth accretion rate. If we 
instead include galaxy mergers (not shown), which increases 
the accretion rate for high-mass galaxies slightly, then the 
results are in fact fully converged for the 50 and 25 h~ 1 Mpc 
box. 

The bottom panel shows that the convergence is better 
at z = 0. 

Although the results of these convergence tests are en- 
couraging, we caution the reader that we cannot exclude the 
possibility that higher resolution simulations would show 
larger differences. For example, iDalla Vecchia &; Schavd 
(2008) found that higher resolution is required to obtained 
converged predictions for the galactic winds. On the other 
hand , it is not clear the resolu tion requirements inferred 
bv lDalla Vecchia fe Schavel (2008) carry over to the present 
simulations, because they used idealised, isolated disk galax- 
ies haloes that started with thinner stellar disks than are 
formed in our simulations and their gas disks were initially 
embedded in a vacuum. 



5.1 Accretion on to haloes 

Figure [3] shows the average specific gas accretion rate on 
to haloes as a function of halo mass for z — 5 (black, top 
curve) to z = (green, bottom curve) for run L050N512. 
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Figure 3. Specific gas accretion rates on to haloes for redshifts 
z = 5 (black, top curve) to z = (green, bottom curve) against 
the total halo mass at the redshift of accretion for the simu- 
lation REF-L050N512. We added the highest mass bins from 
REF_L100N512 to extend the dynamic range. The solid curves 
are for smooth accretion, whereas the dotted curves show the 
specific accretion rate due to mergers with mass ratios greater 
than 1:10. Each mass bin contains at least 10 haloes. Dashed 
lines show fits to da rk matter accretion rates (jPekel et al. 2009b; 
iNeistein et al.l2 006). The specific smooth accretion rate decreases 
with redshift and increases only mildly with halo mass. Except for 
clusters (Mh a i Q ^ 1O 14 M0) at z = 0, halo growth is dominated 
by smooth accretion. 



We also show the high-mass bins from L100N512, which are 
not sampled by L050N512, because of its smaller volume. 
The solid curves show the rates for smooth accretion, the 
dotted curves for mergers with mass ratios greater than 1:10. 
The specific smooth accretion rate decreases with decreasing 
redshift and increases slowly with halo mass. The accretion 
rate, as opposed to the specific accretion rate, is thus roughly 
proportional to the halo mass. 

Our accretion rates are generally in quantitative agree- 
ment with other simulations ([Ocvirk et al . 2008; K eres et al.l 
2005). They also agree well with a fit to an analytic predic- 
tion for dark matter accretion rate s based on the exte nded 
Press-Schechter formalism given by Dekel et al. (2009b) and 
derived by INeistein et al.l (|2QQ6h . This fit is shown by the 
dashed lines. For low- mass haloes (Mhaio ^ 10 12 M©) at 
^ z ^ 2 the accretion rate is lower than predicted. This 
suppression is due to SN winds, as we will discuss below. 
Without SN feedback, the gas accretion rates follow the dark 
matter accretion rates also for low-mass haloes at low red- 
shift. 

The specific gas accretion rate through mergers is much 
lower than the specific smooth accretion rate, except for 
high-mass haloes (Mhaio > 10 14 M©) at z — 0. The domi- 
nance of smooth accretion is consistent with recent work on 
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dark matter accretion, which has shown that mergers with 
mass ratios greater than 1:10 contribute less than 20-30% to 
the total halo accret ion and at least 10-40% of the accretion 
is genuinely diffuse (|Fakhou ri & Ma 2010: lAngulo &; Whitel 
l20ld : iGenel et al.1 l20ld : IWang et alJl20ld ) . 

The power of the OWLS project lies in the fact that 
many simulations have been performed with different pre- 
scriptions for physical processes such as cooling and feed- 
back. They have been described in Section [2] Figure 3] shows 
how the different physical processes affect the specific halo 
accretion rates at z = 2 (black curves) and at z = 
(red curves). At z = 2, we include results from both the 
L025N512 and L100N512 versions of each model in order 
to extend the range of halo masses. These simulations have 
different box sizes and differ by a factor of 64 in mass reso- 
lution. For some models the run with lower resolution and 
larger box size, which is used for the high halo masses, is 
not fully converged, resulting in discontinuities. Only the 
L100N512 simulations were run down to z = 0. 

The first thing to notice is the fact that the halo smooth 
accretion rates are very similar for all models. For the sim- 
ulations with strong winds ( DBLIMFCONTSF VI 6 18) or 
AGN feedback (AGN) the rates are smaller than those for 
the reference model by up to 0.3-0.4 dex. The other mod- 
els differ even less, except for the run without SN feedback 
( NOSN.NOZCO OL) which predicts 0.6 dex higher accretion 
rates at the lowest halo masses. The differences are similar 
at the two redshifts. 

The left panel of Figure U shows the effect of exclud- 
ing metal-line cooling and SN feedback. The simulations 
without SN feedback (NOSN, which is only available for 
L100N512 at z = 2, and NOSN.NOZCOOL) have nearly 
completely flat specific accretion curves. As is the case for 
dark matter accretion rates, there is only a small increase 
with halo mass (see Figure [3]). Including SN feedback with 
^wind = 600 km/s (REF and NOZCOOL) reduces the ac- 
cretion rates for low-mass haloes by up to 0.6 dex (at 
Mhaio = 10 10 M©). It is possible that the true accretion 
rates are the same, but some gas is pushed out of the halo 
before a snapshot is made and is therefore not counted as 
having accreted. A more likely possibility is, however, that 
galactic winds prevented the gas around low-mass haloes 
from accreting. Hence, more gas remains available for accre- 
tion on to massive galaxies which explains why those have 
somewhat higher accretion rates if SN feedback is included. 
Excluding metal-line cooling has no effect on the halo accre- 
tion rates. This indicates that halo accretion rates are not 
set by cooling. 

SN feedback models with higher wind velocities, but 
using the same amount of energy per unit stellar mass, plot- 
ted in the middle panel of Figure 2] give lower accretion 
rates over the range of masses for which the winds are able 
to eject gas. Higher wind velocities keep the feedback ef- 
ficient up to higher ha lo masses ([Dalla Vecchia &; Schavel 
2008; Sch ave et alJl2O10L Haas et al. in preparation). 

The right panel shows the results for models that use 
more energy for feedback than the reference model. A top- 
heavy IMF in starbursts ( DBLIMFCONTSF VI 6 18) reduces 
the accretion rates by up to 0.3 dex. AGN feedback (AGN) 
reduces the accretion rates by up to 0.4 dex for Mhaio > 
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Figure 6. PDF of the lookback time at which the maximum 
past temperature (evaluated at z = 0) was reached for gas ac- 
creted just before z = 5, 4, 3, 2, 1, and in black, blue, purple, 
red, orange, and green, respectively, combining REF-L050N512 
and REF_L100N512. The top a?-axis indicates the corresponding 
redshift. The maximum temperature reached by a gas particle is 
associated with the accretion event, as evidenced by the fact that 
the PDFs peak at the accretion redshift. 



change the halo accretion rates, but very efficient feedback 
from stars and/or AGN can suppress the accretion rates by 
factors of a few. 



5.2 Accretion on to galaxies 

The three lower, thin curves in the top and bottom pan- 
els of Figure [2] show specific smooth accretion rates on to 
galaxies (i.e. the ISM) against total halo mass at z — 2 
and 0, respectively. The specific accretion rate peaks at 
Mhaio ~ 10 12 M©. At z — the peak occurs at slightly 
higher halo masses. For both redshifts the peak falls at 
T v i r rsj 10 6 K, close to the bu mp in the cooling curve due 
to iron (|Wiersma et al.ll2009aT ). At higher halo masses, and 
thus higher virial temperatures, cooling times become long, 
preventing shocked gas from condensing on to galaxies. At 
lower halo masses, feedback prevents gas from entering the 
ISM. For Mhaio ~ 10 12 M the galaxy accretion rate is about 
a factor of two lower than the halo accretion rate, but the 
difference is much larger for other halo masses. 

Figure [5] shows the same as Figure |4] but for accre- 
tion on to galaxies. For models WDENS, DBLIMFCON- 
TSFV1618, and AGN the convergence is poor for the low- 
resolution runs, as can be seen from the discontinuities. Be- 
cause the feedback in these models depen ds on the gas den - 
sity, they are sensitive to the resolution (| Schave et al ll201oh . 
Note, however, that the difference in mass resolution is enor- 
mous, a factor of 64, and that the high-resolution model may 
thus be much closer to convergence than the comparison sug- 
gests. Nevertheless, we caution the reader that the accretion 
rates may be different for higher resolution simulations. 

The variations in the feedback prescriptions result in 
similar differences as for accretion on to haloes, although 
they are generally much larger. Excluding metal-line cooling 
gives slightly different results than halo accretion. 

For Mhaio < 10 13 M© less gas reaches the star formation 
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Figure 4. Specific gas smooth accretion rates on to haloes against total halo mass at z = 2 (top, black curves) and at z = (bottom, red 
curves) for different simulations. The normalization gives an estimate of the time it would take a halo to grow to its current mass at the 
current accretion rate. The curves at low (high) halo masses are for simulations in a 25 (100) /i^Mpc box using 2x512 3 particles. The 
solid curves use a simulation with the reference parameters (REF) and are repeated in all panels. The different simulations are described 
in Section 12.11 Each mass bin contains at least 10 haloes. Each halo contains at least 1000 dark matter particles for the z = 2 curves. 
This is higher than our resolution limit for accretion on to haloes, but it removes the overlap between simulations of different resolution. 
Each halo at z = contains at least 250 dark matter particles. Left panel: Turning off feedback from SNe results in up to 0.6 dex higher 
accretion rates, for low-mass haloes. Middle panel: Increasing the wind velocity causes the specific accretion rate to decrease by up to 
0.2 dex, over the mass range where the feedback is efficient. Right panel: Efficient feedback from a top-heavy IMF or from AGN even 
reduces the accretion rates for the highest halo masses, although the differences between these models are still small, at most 0.3 dex. 
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Figure 5. Specific gas smooth accretion rates on to galaxies against halo mass for different simulations. The curves at low (high) halo 
masses are derived from simulations in a 25 (100) h~ 1 Mpc box using 2 x 512 3 particles. The line styles, colours, and normalization 
are identical to those used in Figure 0] The normalization does not give an estimate for the time it would take the galaxy to grow, 
because we divide by the halo mass, not the galaxy mass. Each mass bin contains at least 10 haloes. In the absence of SN feedback the 
specific accretion rate on to galaxies declines with halo mass, indicating that gas accretes less efficiently on to galaxies in higher-mass 
haloes. Leaving out metal-line cooling decreases the accretion rate most strongly for Mhaio ~ 10 12 Mq. Efficient SN feedback reduces the 
accretion rates substantially for galaxies in low-mass haloes, resulting in a peak in the specific gas accretion rate at Mhaio ~ 10 12 Mq. 
The effects of feedback and metal-line cooling are much stronger for accretion on to galaxies than for accretion on to haloes and can 
result in differences of an order of magnitude. 



threshold (i.e. accretes on to galaxies) without metal-line 
cooling (NOZCOOL), particularly for M ha i ~ 10 12 M at 
z = 2 and M ha io ~ 10 12 4 M at z = 2. The reduction is 
less than 0.4 dex at z = 2, but a full order of magnitude at 
z = 0. For Mhaio > 10 13 M more gas accretes on to galaxies 
at z = if metal- line cooling is excluded. Because there is 
less cooling, less gas accretes on to low- mass galaxies and at 
high redshift. Therefore, there is more gas left to accrete on 
to high mass galaxies at low redshift. 

For simulations without SN feedback the specific galaxy 
accretion rate peaks at Mhaio ~ 10 10 M©, which is two or- 
ders of magnitude lower than when SN feedback is included. 



Without SN feedback, the galaxy accretion rates are a bit 
lower for high halo masses. This could be because much of 
the gas that accreted on to a halo at higher redshift has 
in that case already been accreted on to the galaxy's pro- 
genitors. Alternatively, the galaxy accretion rates could be 
higher in the presence of SN feedback due to the increased 
importance of recycling. If winds are able to blow gas out 
of the galaxy, but not out of the halo, as may be the case 
for high-mass haloes, then the same gas el ements may be ac- 
crete d on to the galaxy more than once ([O ppenheimer et "all 
[20l3). For M ha i < 10 11 M at z = 2 and M ha io < 10 12 M 
at z = 0, on the other hand, the accretion rates are higher, 
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because more gas accretes on to the halo and because there 
is no feedback to stop halo gas from accreting on to the 
galaxy. 

Galaxy accretion rates are thus much more sensitive to 
metal-line cooling and to SN and AGN feedback than halo 
accretion rates. The difference between models can be as 
large as 1 dex. Feedback processes determine the halo mass 
for which galaxy formation is most efficient. 



6 HOT AND COLD ACCRETION ON TO 
HALOES 

Figure [6] demonstrates that most of the gas reaches its max- 
imum temperature around the time it is accreted on to a 
halo. Here the probability density function (PDF) of the 
lookback time at which the gas reaches its maximum tem- 
perature, evaluated at z — 0, is shown for baryonic particles 
that were accreted as gas particles on to haloes at different 
redshifts. The vertical dotted lines show the times at which 
the gas was accreted on to a halo. The fact that the PDFs 
peak around the accretion redshifts shows that the maxi- 
mum temperature is usually related to the accretion event. 

Some of the gas reaches its maximum temperature sig- 
nificantly later than the redshift at which it was accreted. 
When two galaxies merge, gas can shock to higher temper- 
atures. Gas can also be affected by winds resulting from 
SN feedback. Because we are primarily interested in gas ac- 
cretion, we will from now on evaluate the maximum past 
temperature of the gas at the first available output redshift 
after accretion. 

Gas can be accreted cold on to haloes with well devel- 
oped virial shocks if its density is sufficiently high, as can 
for example be the case in filaments. To illustrate this we 
show in the left two columns of Figure [3 the gas overdensity 
and the temperature in a cubic region of 1 h~ x Mpc (comov- 
ing) centred on three example haloes at z — 2 taken from 
the full sample of 12768 haloes in the high-resolution ref- 
erence simulation (REF_L025N512). The haloes have total 
masses M ha io « 10 12 5 , 10 12 , and 10 115 M (from top to bot- 
tom), corresponding to comoving virial radii of 379, 264, and 
170 /i -1 kpc, respectively, shown as white circles in the tem- 
perature plots. Their virial temperatures are T v i r « 10 6-3 , 
10 6 , and 10 5 ' 7 K, respectively. The colour scales are the 
same for all three haloes. We can immediately see that 
the average temperature increases with halo mass. Hot gas, 
heated either by accretion shocks or SN feedback, extends 
to several virial radii. Without SN feedback, the hot gas 
would trace the virial radius more accurately for the two 
lowest mass haloes, as can be seen for the 10 12 Mq halo in 
Figure fT3l 

Most of the gas around the 10 12-5 Mq halo has been 
heated to temperatures above 10 6 K and the halo will get 
most of its gas through hot accretion. Cold streams do pen- 
etrate the virial radius, but they seem to break up as they 
get close to the centre. Even so, a number of small, dense 
clumps do survive and remain relatively cold. Such cold 
clum ps originating from filame ntary gas were als o stu died 
bv lSommer-Larsenl (|2006h and lKeres fe Hernauistl (|2009h . 

The 10 12 M© halo is located at the intersection of three 
filaments. The gas in the filaments is denser and colder 
than in the surrounding medium. The cold streams become 



narrower as they get closer to the centre. They are com- 
pressed by the high pressure, shock-heated gas around them 
(|Keres et al.l l2009ah . Cold streams bring gas directly and 
efficiently to the inner halo. 

The 10 115 M halo is embedded in a single filament. 
It has the lowest virial temperature, so the hot gas is much 
colder than in the highest mass halo. Cold streams are most 
prominent and broadest in this halo. This halo will get most 
of its gas through cold accretion. 

The right two columns of Figure shows zooms of the 
central 250 h' 1 kpc (comoving). The 10 12 5 M halo con- 
tains a large number of cold clumps. The galaxies in the 10 12 
and 10 11 ' 5 M haloes are being fed by cold, dense streams. 
All of these galaxies have formed discs. The galaxy in the 
10 12 ' 5 Mq halo has clear spiral arms and a bar- like struc- 
ture. The galaxy in the 10 12 M© halo has a very small disc 
with cold gas around it, which looks more disturbed. The 
galaxy in the 10 11-5 Mq halo is fairly large with a lot of cold 
material accreting on to it. 

The maximum temperature reached by shock-heated 
gas is expected to scale with the virial temperature of the 
halo. However, we do not expect the ratio of T max and T v i r 
to be exactly unity, because of departures from spherical 
symmetry, adiabatic compression after virialization, and the 
factor of a few difference between different definitions of T v i r . 

Dividing T max by T v i r would take out the redshift and 
halo mass dependence of the virial temperature. We calcu- 
late the virial temperature as follows 

_ ( G 2 Hln m l^ \ 1/3 M m H 2/3 , , 
Jv,r " 1, 54 ) k B M ^ L + Z >' 

where G is the gravitational constant, Ho the Hubble con- 
stant, fi the mean molecular weight, mu the mass of a hy- 
drogen atom, and ku Boltzmann's const ant0. 

6.1 Dependence on halo mass 

Figure [8] illustrates the dependence of the maximum past 
temperature on halo mass. Shown are scatter plots of the 
maximum past temperature reached by gas accreting on to 
haloes just before z — 2 (top panel) and z — (bottom 
panel) against the mass of the halo at these redshifts. The 
logarithm of the accreted gas mass per (/i _1 Mpc) 3 in a pixel 
is used for colour coding. The virial temperature is indicated 
by the black line. 

A clear bimodality is visible for accretion at z = 2, with 
a minimum at T max ~ 10 5-5 K, indicated by the dotted line. 
This minimum coincides wit h a maximum at T 10 5 ~ 5-5 K 
in the cooling function (e.g. IWiersma et a l. 2009a). For gas 
accreted in the hot mode, which includes most of the gas 
accreting on to high-mass haloes, the maximum past tem- 
perature is within a factor « 3 of the virial temperature and 
displays the same dependence on mass. The temperature of 
the gas accreted in the cold mode is independent of halo 
mass. 

3 This definition is a factor 2/3 lower than t he virial tem perature 
used by some other authors (e.g. Barkana &; Loeb l200lh . 
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Figure 7. Gas overdensity (first and third columns) and temperature (second and fourth columns) in a cubic region of 1 h~ x comoving 
Mpc (first and second columns) and 250 h~ x comoving kpc (third and fourth columns) centred on haloes of M^aio ~ 10 12,5 , 10 12 , and 
10 11 - 5 Mq (from top to bottom) at z = 2 fo r simulation REF_ L025N512. The white circles indicate the virial radii of the haloes, as 
computed using the overdensity criterion from lBrvan fc Normanl |l998). Cold, dense streams bring gas to the centre. The temperature of 
the hot gas increases with halo mass. Hot accretion dominates for high-mass haloes, cold accretion for low-mass haloes. The galaxies in the 
centres of these haloes are discs, surrounded by cold gas. This cold gas is in clumps (Mh a io ~ 10 12-5 ), disrupted streams (Mh a i ~ 10 12 ), 
or smooth streams (Mh a i Q « 10 11-5 ). 



At z = 0, the lowest T ma x values are higher than at 
z = 2. This shift occurs because the density, and hence the 
cooling rate, increases with redshift. At fixed halo mass, the 
highest Tmax values are lower at z = 0, because the virial 
temperature of a halo at fixed mass decreases with decreas- 
ing redshift, as can be seen from Equation [TJ For haloes with 
virial temperatures T v i r < 10 5 K it becomes impossible to 
tell from this plot whether or not the gas has gone through 
a virial shock because the virial temperature is similar to 
the maximum past temperature reached by gas accreting in 
the cold mode. This makes it difficult to separate hot and 
cold accretion for low-mass haloes at low redshift. As we 
will show below, separating hot and cold accretion using a 
fixed value of T m ax/T v ir is more difficult than using a fixed 
value of T max , because the minimum in the distribution is 
less pronounced and because it evolves (|Keres et al.ll2005r ). 
In most of this paper, we will therefore use a fixed maximum 
temperature threshold of T ma x = 10 5-5 K. 



The relativ e importance of hot accretion increases with 
halo mass (e.g. lOcvirk et al J 12003 ). The top panel of Fig- 
ure [9] shows the fraction of gas smoothly accreting on to 
haloes in the hot mode, just before z = 2 for simulations 
REF.L100N512, REF.L050N512, and REF.L025N512. The 
bottom panel shows this for accretion just before z — 
for simulations REF.L100N256, REF.L100N512, and 
REF-L050N512. A particle accreted just before z = 2 is con- 
sidered to have been accreted hot if T ma x(z = 2) ^ 10 5-5 K. 
The error bars show the la halo to halo scatter. 

We have checked, but do not show, that the results are 
fully converged with box size for fixed resolution. In each 
panel the three simulations span a factor 64 in mass reso- 
lution. The hot fraction decreases slightly with increasing 
resolution, but the differences are very small. This slight 
decrease could arise because higher density regions inside 
clumps and filaments are better sampled with increasing res- 
olution, leading to higher cooling rates in cold gas. 
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Figure 8. Maximum past temperature at z = 2 (top panel) 
and z = (bottom panel) against total halo mass of the 
gas smoothly accreted on to haloes for models REF_L025N512 
and REF_L100N512 (top panel) and REF.L050N512 and 
REF.L100N512 (bottom panel). The logarithm of the total gas 
mass per (/i -1 Mpc) 3 in a pixel is used for colour coding. The 
solid line indicates the virial temperature of the halo. The dot- 
ted line shows T m ax = 10 5-5 K, where there is a minimum in the 
mass per pixel at z = 2. The temperature of gas that is accreted 
hot scales with the virial temperature. For low-mass haloes, the 
temperatures of hot and cold accreted gas are comparable. Hot 
mode accretion is more important for higher halo masses and 
lower redshifts. 



Hot mode accretion becomes indeed more important for 
higher mass systems. This is expected because only suffi- 
ciently massive haloes are capable of providing pressure sup - 
port for a stable virial shock (e.g. iBirnboim &; DekellfeoOcfl ). 
The median hot fraction (not shown) behaves similarly to 
the mean, although it is smaller for low-mass haloes. 

There is no sharp transition from cold to hot accretion. 
At z = 2, the hot mode accretion increases from 20% to 80% 
when halo mass increases from 10 11 to 10 13 M©. The proper- 
ties of galaxies are therefore not expected to change suddenly 
at a particular halo mass as has been assumed in some semi- 
analy tic models (e.g. ICattaneo et al . 2008; Cr oton fc Farrarl 
120081 ). 

In Figure [10] we compare this result to that obtained 
when we define hot mode accretion using a maximum tem- 
perature threshold that depends on the virial temperature. 
Using T max = Tvir shows that, even for the most massive 
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Figure 9. Average fraction of the gas, smoothly accreted on 
to haloes between z = 2.25 and z = 2 (top panel) or between 
z = 0.125 and z = (bottom panel), that has maximum past 
temperature T max ^ 10 5,5 K. The different curves are from simu- 
lations of the same reference model but spanning a factor of 64 in 
mass resolution for each panel. The error bars show the la halo 
to halo scatter for simulation REF_L100N512. Each mass bin 
contains at least 10 haloes. Arrows correspond to the adopted 
resolution limit for accretion on to haloes. Cold mode accretion 
dominates for M^aio < 10 12 Mq, but the transition is very grad- 
ual. 



haloes, only about 40% of the gas reaches a temperature 
higher than T v i r at z — 2. At z — about 70% of the gas ac- 
creting on to massive haloes reaches T v i r . Because gas going 
through a virial shock may only heat to a factor of a few be- 
low Tvir, this definition does not discriminate well between 
hot and cold accretion. If we decrease the critical T ma x/Tvir, 
then the hot fraction gets close to our previously determined 
value for the high-mass haloes. For low-mass haloes, how- 
ever, this results in a sharp upturn of /hot as it must ap- 
proach unity if the threshold maximum temperature falls 
much below 10 5 K. 

The hot fraction therefore depends very much on the 
definition of the temperature threshold. Depending on the 
definition, gas accreted cold may in fact have experienced a 
virial shock. For haloes with T v i r ^> 10 5 K, however, we can 
safely separate hot and cold accretion and trust the result 
that a larger fraction of the gas goes through a virial shock 
for higher- mass haloes. 
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Figure 11. Hot fraction for accretion on to haloes just before 
z = 5, 4, 3, 2, 1, and against halo mass at the same red- 
shift. The curves at low halo masses are obtained from simulation 
REF-L050N512. At the high-mass end we have added curves for 
simulation REF_L100N512 to extend the dynamic range. Gas is 
considered to have been accreted hot if it has T max ^ 10 5-5 . The 
dotted curves show the hot fraction including both smooth accre- 
tion and mergers. Each mass bin contains at least 10 haloes. For 
a fixed halo mass, hot accretion tends to be more important at 
lower redshift. 



differences are negligibly small. Indeed, the hot fraction for 
gas accreted in mergers (not shown) is nearly the same as 
that for gas accreted smoothly. 



Figure 10. The average fraction of the gas, accreted on to haloes 
just before z = 2 (top panel) or z = (bottom panel), that has 
maximum past temperature above a certain temperature thresh- 
old. The results are for simulation REF_L050N512 with the high- 
est mass haloes from REF.L100N512. The thresholds for the 
blue, purple, red, orange, and green curves (bottom to top) are 
T m ax/T v i r = 1, lO" 25 , lO" 05 , 1CT - 75 , and 0.1, respectively. 
The black curve shows the average fraction with T ma x ^ 10 5-5 K 
and is identical to the black curve in Figure [9] Each mass bin 
contains at least 10 haloes. The hot fraction depends strongly on 
the choice we make for the threshold, particularly for lower-mass 
haloes. 

In the rest of this paper, we will define hot mode accre- 
tion using a fixed maximum past temperature threshold of 
T max = 10 5 - 5 K. 



6.2 Smooth accretion versus mergers 

So far we have only looked at 'smooth' accretion, i.e. we ex- 
cluded mergers with a mass ratio greater than 1:10. We could 
have chosen not to exclude mergers, because the gas reser- 
voir of a halo can also grow through mergers. Even though 
mergers with ratios greater than 1:10 contribute only ~ 10% 
of the total gas accretion for Mhalo < 10 14 Mq (see Figure[3]), 
it is interesting to investigate the differences for hot and cold 
accretion. 

Figure [TT] compares the hot fraction of smooth accre- 
tion (solid curves) and all accretion (dotted curves), which 
takes into account both smooth accretion and mergers. The 



6.3 Dependence on redshift 

The fraction of gas that is accreted on to haloes in the hot 
mode also depends on redshift, as can be seen from Fig- 
ure [TT] where the hot fraction is plotted against halo mass 
for different redshifts. For a given halo mass, the hot frac- 
tion increases with time between z — 5 and z — 1. Below 
z — 1 the rate of evolution slows down, presumably because 
structure formation slows down, and the sign of the evolu- 
tion may reverse for low-mass haloes due to the decrease of 
Tvir with time. 

At high redshift the proper density of the Universe is 
higher and the cooling time is therefore shorter (t coo i oc 
p~ x oc (1 + z)~ 3 ). The Hubble time is also shorter, but its 
dependence on redshift is weaker (tn oc H~ x oc (Q m (l + 
z) 3 + n A )" 1/2 , so tn oc (1 + z)~ 3/2 at high redshift and t H 
is independent of redshift at low redshift). The dynamical 
time at fixed overdensity has the same redshift dependence 
as tn- Hence, cooling is m ore efficient at higher redshifts (e.g. 
iBirnboim fc Dekelll20o3 ) and the hot fraction will thus be 
lower for a fixed virial temperature. For a fixed halo mass, 
the evolution in the hot fraction is smaller, because T v i r is 
higher at higher redshifts. The environment can also play 
a role. At high redshift (z > 2), 10 12 M© haloes are rare 
and they tend to reside in highly overdense regions at the 
intersections of the filaments that make up the cosmic web. 
Haloes of the same mass are, however, more common at 
low redshift and form in single filaments, with more average 
densities. Hence, cold streams may be able to feed massive 
haloes at hig h redshift, whereas this may not be possibl e at 
low redshift (jKeres et al.ll2005l : iDekel fe Birnboimll2006T ). 
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Figure 12. Specific smooth accretion rates of gas on to 
haloes against total halo mass just before redshifts z = 5 
(black, top curve) to z = (green, bottom curve) for the 
simulation REF-L050N512, with the highest mass bins from 
REF_L100N512. The solid and dashed curves are the rates for 
hot (T max > 10 5 - 5 K) and cold (T max < 10 55 K) accretion, re- 
spectively. Each mass bin contains at least 10 haloes. The specific 
hot accretion rate increases with halo mass, whereas the specific 
cold accretion rate decreases for Mh a i Q > 10 12 Mq. 



The specific accretion rate of gas on to haloes increases 
mildly with halo mass for < z < 5, see Figure [3] Splitting 
this into hot and cold specific accretion rates reveals a steep 
increase with halo mass for hot accretion, as shown in Fig- 
ure U21 The specific cold accretion rate decreases with halo 
mass for M halo > 10 12 M . 

From observations we know that the specific star forma- 
tion rate, i.e. the SFR divided by the st ellar mass, declines 



with both time and ste ll ar mass (e.g. iBrinchmann et al 
2004; iBauer et al.1 120051 : iFeulner et all 120051 : IChen et al 
2009). The decline in the specific SFR may be related to 



the decline in the specific cold accretion rate. However, the 
decline is much stronger in the observations and present for 
the entire stellar mass range probed. 

6.4 Effects of physical processes 

As discussed in Section [5l feedback from SNe and AGN de- 
creases the halo accretion rate, while metal-line cooling has 
very little effect. In this Section we will discuss the influ- 
ence of cooling and feedback on the hot and cold accretion 
fractions at z — 2 and z = 0. 

To investigate the influence of galactic winds driven by 
SN feedback, we ran simulations with no SN feedback at all 
and with more effective galactic wind models, that can eject 
gas from more massive haloes. In other simulations, cooling 
rates are computed assuming primordial abundances. We 
have also included AGN feedback in one simulation. 

To illustrate the effect these different processes have 
on an individual halo, Figure [13] shows the same 10 12 M 
halo as was shown in the middle panel of Figure for five 
different simulations. From top to bottom are shown: no 
SN feedback and no metal- line cooling (NOSN.NOZCOOL); 
no metal- line cooling (NOZCOOL); reference SN feedback 
(REF); density dependent SN feedback (WDENS); refer- 
ence SN feedback and AGN feedback (AGN). The colour 
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Figure 13. Gas overdensity (left) and temperature (right) in a 
1 comoving Mpc box centred on a halo of 10 12 Mq at z = 2. 
We show the same halo as was shown in the middle panel of Fig- 
ure[71 but now for five different models. The white circles indicate 
the virial radius of the halo. From top to bottom: no SN feed- 
back and no metal-line cooling (NOSNNOZCOOL), no metal- 
line cooling (NOZCOOL), reference SN feedback (REF), density 
dependent SN feedback (WDENS), reference SN feedback, and 
AGN feedback (AGN). These simulations used 2 x 512 3 particles 
in a 25 h~ 1 Mpc box. The structure of the cold streams changes, 
but they exist in all simulations. The hot gas extends to larger 
radii and has a higher temperature if feedback is more efficient. 
The structure inside the halo is clearly different in different sim- 
ulations, but cold gas is always present outside the disc. 
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coding shows gas overdensity in the left panels and temper- 
ature in the right panels. The stronger the feedback (top 
row: weakest feedback, bottom row: strongest feedback) the 
more fragmented the streams become. All models predict the 
presence of cold, dense gas throughout much of the halo. The 
diffuse halo gas is heated to higher temperatures in simula- 
tions with strong SN or AGN feedback than in the reference 
simulation. The radius out to which gas is heated increases 
for strong feedback models. 

Figure [14] shows how the fraction of the gas that ac- 
cretes smoothly on to haloes just before z — 2 (top pan- 
els) and z = (bottom panels) in the hot mode, i.e. 
with T max ^ 10 5 ' 5 K, depends on the physical processes 
that are modelled. Even though the hot fractions are not 
completely converged at the low-mass end for some of the 
L100N512 runs, the trend with halo mass and the effect 
of the variations in the simulations are robust. The haloes 
above M ha i « 10 12 M in the L025N512 runs (not plot- 
ted) show an increase in the hot fraction as steep as the 
L100N512 runs. Before discussing the differences between 
the models, we stress that these differences are small. The 
fraction of the gas that accretes on to a halo of a given mass 
and at a given redshift in the hot mode can thus be robustly 
estimated. It is insensitive to uncertainties in the baryonic 
physics, such as radiative cooling and feedback from star 
formation and AGN. 

We first focus on the left panels, which compare sim- 
ulations with and without metal-line cooling and with and 
without SN feedback. For high-mass haloes the fraction of 
gas accreting in the hot mode is a little bit higher without 
metal-line cooling, because the gas will reach higher temper- 
atures if the cooling rates are lower. The effect is, however, 
small. 

The hot fraction depends only slightly on the specific 
feedback model used if the energy per unit stellar mass is 
kept fixed, as can be seen from the middle panels of Fig- 
ure 1141 The hot fraction is a bit lower for the model with 
the most effective feedback (WDENS). 

The right panels show that AGN feedback decreases the 
hot fraction, though the effect is not large and limited to 
high-mass (^> 10 11 M©) haloes. More effective stellar feed- 
back as a result of a top-heavy IMF in starbursts (DBLIMF- 
CONTSFV1618) reduces the hot fraction more for low-mass 
< 10 12 M© haloes. These results suggest that the hot accre- 
tion mode is slightly more affected by feedback than the cold 
mode. Because the hot gas is less dense than the cold gas 
and because it spans a greater fraction of the sky as seen by 
the galaxy, it is more likely to be affected by feedback. It has 
been argued by several authors that AG N feedback would 
therefore couple mostly to the hot gas (jKeres et alJ 120051 : 
iDekel fc Birnboiml l 2006 V Our results show that this effect 
is small, although we will show below that it is important 
for accretion on to galaxies residing in high-mass haloes at 
low redshift. The main conclusion is that the fraction of the 
gas that accretes on to haloes in the hot mode is insensitive 
to feedback from SNe and AGN. 



7 HOT AND COLD ACCRETION ON TO 
GALAXIES 

A significant fraction of the gas that accretes on to haloes 
may remain in the hot, low-density halo. This gas never cools 
and it will not get into the inner galaxy and contribute to 
the SFR. This diffuse gas may be easily pushed out of the 
halo by galactic winds. It is therefore of interest to look not 
only at the gas that accretes on to haloes, but also at the 
gas that actually accretes on to galaxies. After all, it is only 
the gas that is accreted on to galaxies that is available for 
star formation. 

Gas is by definition cold when it accretes on to a galaxy, 
because it must cool below T — 10 5 K to be able to join 
the ISM. However, by Using Tmax we are probing the entire 
temperature history of the gas and not just the temperature 
at the time of accretion. Gas accreting on to a galaxy in the 
hot mode has been hot in the past (usually when it accreted 
on to the halo), but was able to cool down and reach the 
central galaxy. Gas accreting on to a galaxy in the cold mode 
has never been hot in the past. 

In Figure [15] we show the hot fraction for accretion on 
to haloes (dotted curve), accretion on to the ISM (solid 
curve), and for stars formed (dashed curve) just before z — 2 
(top panel) and just before z = (bottom panel) for the 
50 /i -1 Mpc reference simulation. To illustrate the conver- 
gence, we also show the hot fraction for accretion on to the 
ISM for simulations with different resolutions. We show all 
results down to halo masses corresponding to 100 dark mat- 
ter particles, i.e. 1 dex below our resolution limit. 

For high-mass haloes (Mhaio > 10 12 M©) hot mode ac- 
cretion is less important for accretion on to the ISM, and 
therefore for star formation, than for accretion on to haloes. 
For T vir > 10 6 K the hot fraction remains approximately 
constant with mass, even though the hot fraction for halo 
accretion becomes larger. The lower hot fractions arise be- 
cause the temperature of the hot gas increases with the 
virial temperature and for th ese temperatures hot ter gas 
has a longer cooling time (e.g. IWiersma et al.ll2009al ). mak- 
ing it less likely to enter the galaxy. In reality, cold, dense 
clouds may be disrupted more easily than in SPH simula- 
tions , which could push the hot fraction up ([Agertz et al.1 
M). 

On the other hand, for low-mass haloes (T v i r < 10 5-5 K) 
the hot fraction for accretion on to galaxies is higher than 
for accretion on to haloes. The virial temperature of these 
haloes is so low that the maximum temperature will only be 
above 10 55 Kifthe gas was heated by SN feedback. This can 
happen after the gas has accreted on to the halo, but before 
the gas joins the ISM, explainin g the higher hot fract i on for 
accretion on to galaxies. Indeed. lO ppenheimer et all ([2010) 
have shown that the re-accretion of gas that has been ejected 
by galactic winds can be important. The hot fraction for 
low-mass haloes is lower in simulations without feedback, 
confirming our interpretation (see Figure [17}. 

As expected, for recently formed stars the hot fraction 
is comparable to that for gas that recently accreted on to the 
ISM, although it is a bit lower. It is slightly lower because 
of the time delay between accretion and star formation. The 
gas from which the stars were formed was typically accreted 
at higher redshift and on to lower-mass haloes, which corre- 
sponds to lower hot fractions. 
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Figure 14. The average fraction of the gas that accreted smoothly on to haloes just before z = 2 (top panels) and z = (bottom 
panels), that reached maximum past temperatures above 10 5 ' 5 K (by z = 2 and z = 0, respectively) is plotted as a function of total halo 
mass. The line styles are identical to those used in Figure [4] The curves at the high-mass end are from simulations L100N512. At z = 2 
the curves at the low-mass end are from simulations L025N512, resulting in a factor 64 higher mass resolution. Each mass bin contains 
at least 10 haloes. At z = 2, each halo contains at least 1000 dark matter particles. This is higher than our resolution limit for accretion 
on to haloes, but it removes the overlap between simulations of different resolution. The differences between the simulations are small. 
Efficient feedback generally reduces the hot fraction, indicating that hot gas is more vulnerable to feedback than cold gas. 



Our lowest resolution simulations underestimate the hot 
fraction somewhat at their low-mass ends, but this regime 
is excluded by our adopted resolution limit, as indicated by 
the arrows. However, for z = the results are also not fully 
converged at the high-mass end (Mhaio > 10 12 M©). For 
high-mass haloes the hot fraction decreases with increas- 
ing resolution, as we also saw for accretion on to haloes 
(see Figure [9}, although the effect is much larger for accre- 
tion on to galaxies. This is likely because higher densities 
can be reached with higher resolution. Cold mode accretion 
may thus be somewhat more important for fuelling massive 
galaxies than our simulations suggest. 

The evolution of the hot fraction of gas accreted on to 
the ISM is shown in Figure [16] from z = 5 (bottom, black 
curve) to z = (top, green curve) as a function of halo mass. 
The trend with redshift is the same as for accretion on to 
haloes (see Figure [TT]) . The hot fraction increases with de- 
creasing redshift for all halo masses, except when the virial 
temperatures fall below 10 5 ' 5 K, the value we use to sepa- 
rate the hot and cold accretion modes, which happens for 
M ha i < 10 11,5 M at z < 1. The hot fraction for accretion 
on to galaxies increases less steeply with halo mass than 
for accretion on to haloes. This is due to the fact that the 



cooling time of the hot gas increases for higher-mass haloes, 
making it less likely that hot gas reaches the central galaxy. 

For haloes with T v i r > 10 6 K, the hot fraction of gas 
that reaches the ISM is much lower than the hot fraction 
of all the gas that accretes on to the halo. At high redshift 
[z > 4) cold accretion is the dominant mode for feeding 
galaxies. At lower redshift (z < 2) hot and cold accretion 
are comparable. Except for high- mass haloes at low redshift 
(Mhaio > 10 13 Mq at z = 0), hot mode accretion is less 
important for feeding the central galaxy than cold mode ac- 
cretion. It is always less important for the growth of galaxies 
than it is for the growth of haloes, though it is never negli- 
gible. We found the hot fraction of recently formed stars to 
be slightly lower, because it takes some time to convert the 
gas into stars. Cold mode accretion is most important for 
the total build-up of stellar mass in galaxies. 

We have not included galaxy mergers. Including merg- 
ers preferentially brings in gas accreted in the cold mode, 
because that gas was already part of the ISM. This reduces 
the hot fraction slightly for high-mass haloes at low red- 
shift. Cold mode accretion is therefore also the main mode 
for galaxy growth in this case. 
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Figure 15. The dotted curve shows the average fraction of gas 
smoothly accreted on to haloes at z = 2 (top panel) and z = 
(bottom panel), that has maximum past temperatures above 
T max = lO 5-5 K. The solid curves show the average fraction of the 
gas smoothly accreted on to the ISM with maximum past tem- 
peratures above T max = 10 5,5 K. The dashed curves show the 
fraction of the stars formed in the same redshift intervals from 
gas with maximum past temperatures above T max — lO 5 5 K. 
These reference simulations span a factor of 64 in mass resolution 
in each panel. Each mass bin contains at least 10 haloes. Ar- 
rows correspond to the adopted resolution limit for accretion on 
to galaxies. For high-mass haloes, the hot fractions for accretion 
on to haloes and galaxies diverge. The hot fraction for recently 
formed stars follows the same trend as for accretion on to the 
ISM. 



7.1 Effects of physical processes 

We showed the effect of feedback, metal-line cooling, and 
cosmology on galaxy accretion rates in Section [5] Feedback 
reduces these rates, while including metal-line cooling in- 
creases them. In this Section we discuss their consequences 
for the relative importance of hot and cold accretion. 

Figure [T71 shows the hot fraction of the gas that reaches 
the ISM between z = 2.25 and z = 2 as a function of halo 
mass for the same simulations as were shown for accretion 
on to haloes in Figure [14] Even though the hot fractions 
are not completely converged for the L100N512 runs, the 
trend with halo mass and the effect of the variations in the 
physics are robust. As was the case for accretion on to haloes, 
the trends are the same for all simulations, although the 
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Figure 16. Hot fraction against halo mass for accretion on to 
the ISM just before z = 5, 4, 3, 2, 1, and shown from bot- 
tom to top by the curves in black, blue, purple, orange, red, and 
green, respectively. The curves at low halo masses are obtained 
from model REF_L050N512. At the high-mass end we have added 
curves for REF-L100N512 to extend the dynamic range. Each 
mass bin contains at least 10 haloes. The hot fraction increases 
with redshift, except for haloes with virial temperatures which 
fall below 10 5-5 K, the value we use to separate hot and cold ac- 
cretion (M halo < 10 11,5 Mq and z < 1). The hot fraction varies 
less strongly with halo mass than for accretion on to haloes. 



differences between models are larger for accretion on to the 
ISM than for accretion on to haloes. At z = 2, the hot mode 
is less important for accretion on to the ISM, and therefore 
less important for star formation, than the cold mode. For 
the halo mass range shown at z = (> 10 12 Mq), the hot 
mode is slightly more important than the cold mode. 

We can observe the difference between simulations with 
and without metal-line cooling and with and without SN 
feedback in the left panels. The simulation without metal- 
line cooling (NOZCOOL) has a slightly smaller hot fraction 
than the simulation with metal- line cooling (REF) at z = 2 
and for Mhalo < 10 12-5 Mq at z = 0. Because the cooling 
times are longer in the absence of metal-line cooling, less 
hot gas is able to reach the high densities that define the 
ISM. For accretion on to haloes we found the opposite effect 
(see Figure [14]), because the lower cooling rate increases the 
maximum temperature reached by gas accreted on to haloes 
and hence increases the corresponding hot fraction. For high 
mass haloes (Mhalo > 10 13 Mq) at z = 0, the hot fraction is 
much higher in the absence of metal- line cooling. In Figure [5] 
we have seen that metal-line cooling strongly increases the 
accretion rates on to galaxies in lower mass galaxies. Thus, 
without metal-line cooling, less hot gas accretes on to low- 
mass haloes, leaving more hot gas to accrete on to their 
higher mass descendants. 

We can compare the simulation without metal-line cool- 
ing (NOZCOOL) to the one without SN feedback and with- 
out metal- line cooling (NOSN_NOZCOOL). In the simula- 
tions with feedback a higher fraction of the gas accreted on 
to the ISM has been hot, presumably because some of the 
gas accreted cold on to haloes was heated by outflows be- 
fore it joined the ISM (for the last time). At the high-mass 
end at z = the differences between the curves parallels 
those in the specific accretion rates (Figure HJ). This sug- 
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Figure 17. The average fraction of the gas accreting on to the ISM at z = 2 (top panels) and z = (bottom panels), that has maximum 
past temperatures T ma x ^ 10 5-5 K. The line styles are identical to those used in Figure ITTI The curves at the low and high-mass ends 
are from simulations L025N512 and L100N512, respectively, which differ by a factor 64 in mass resolution. Each mass bin contains at 
least 10 haloes. Differences in results from different simulations are mostly small, except for M^aio > 10 13 Mq at z = if metal-line 
cooling is ignored or AGN feedback is added. 



gests that here SN feedback increases the hot fraction be- 
cause it prevents the accretion of hot gas on to lower mass 
progenitors, leaving more gas available to cool in high- mass 
galaxies, where the SN feedback is inefficient. The effect of 
SN feedback on the hot fraction at z — 2 is much smaller 
if metal- line cooling is included (NOSN). This may indicate 
that the increase in the cooling rates due to the metals car- 
ried by the winds compensates for the extra shock heating. 

The middle panels show the result for simulations with 
different values for the SN feedback parameters, but the 
same amount of SN energy per unit stellar mass formed. 
The hot fractions are similar, though slightly higher if the 
feedback model is more efficient ( WML1 V848 for Mhaio > 
10 11 M and WDENS for M ha i > 10 12 M at z = 2). 

Results for simulations with a prescription for SN feed- 
back that use more energy and for a model including AGN 
feedback are shown in the right panels. Perhaps surprisingly, 
at z = 2 the simulation including AGN feedback (AGN) 
gives results similar to the reference simulation. The accre- 
tion rate on to the galaxy is suppressed by up to 1 dex, as 
can be seen in Figure but the hot fraction is almost the 
same. At z — the hot fraction is somewhat lower, though 
still above 40%. 

In Figure [TJ] we showed that the hot fraction for the 
gas accreting on to haloes does not change much if we vary 
the prescriptions for feedback and radiative cooling. To first 



order, the same conclusion is true for the gas that accretes 
on to the ISM and thus becomes star forming, although the 
differences are larger. At z = 2 the hot fraction decreases 
significantly if both SN feedback and metal-line cooling are 
turned off. For high-mass galaxies (Mhaio > 10 13 M©) at 
z = the importance of the hot mode increases substantially 
if metal-line cooling is excluded and decreases significantly 
if AGN feedback is added. The specific implementation of 
the SN feedback does not have a large effect on the fraction 
of the gas accreting on to galaxies that has a maximum past 
temperature T max ^ 10 5 5 K. 



8 COMPARISON WITH PREVIOUS WORK 

We can co m pare our res ults to the pioneer i ng wo rk of 
iKeres et al.1 (|2005l . l2009sJ ). and lOcvirk et al] <|2008h . All 
these studies used 250,000 K as the critical temperature to 
separate hot and cold accretion. We used 10 5 ' 5 = 316, 228 K, 
but our results would have been very similar if we had used 
250 000 K. 

lOcvirk et al.1 (|2008h use an adaptive mesh refinement 
(AMR) code and include metal enrichment and cooling, but 
only weak SN feedback and no AGN feedback. They use 
only a single simulation, whose resolution is similar to our 
L025N512 runs. Their simulation was only run to z ~ 1.5, so 



© 2010 RAS, MNRAS 000, [HS] 



20 F. van de Voort et al 



we can only compare to our high-redshift results. Their box 
size of 50 /i -1 Mpc is too small to sample haloes with mass > 
10 13 M©. As their simulation is not Lagrangian, they cannot 
trace the gas back in time, and were forced to separate the 
cold and hot modes based on current temperatures. 

For z — 2 — 3 they predict that the hot fraction for 
accretion on to haloes reaches 0.5 at Mhaio ~ 10 11 ' 5 M©, 
whereas we find Mhaio ~ 10 12 M©. For higher redshifts we 
predict somewhat higher hot fractions. We predict a gradual 
increase in the transition mass with redshift, but see no sign 
of the sudden change between z — 3 and 4 that they found. 

For ac cretion on to galaxie s, which they measure at 0.2 
virial radii. lOcvirk et al.l (|2008l ) also do not predict a sudden 
change with redshift. As the gas at this radius is not neces- 
sarily star forming and has not necessarily cooled down to 
T ^ 10 5 K, it is not possible to make a completely fair com- 
parison. They find no change with redshift, whereas we find 
strong evolution. The hot fraction increases up to their high- 
est halo mass, reaching values close to unity for 10 12 ' 5 M©. 
In our simulations, however, the hot fraction for accretion 
on to galaxies is approximately constant and never exceeds 
0.4 at z ^ 2. As we are using the SPH technique, it is possi- 
ble that our simulations suffer from in-shock cooling, which 
would lead to an underestimate of the hot fraction. How- 
ever, we find that the hot fraction in fact decreases with 
increasing resolution, the opposite of what we would expect 
if this were an important effect. Perhaps the difference is 
due to the fact that we do not update the maximum past 
temperature once the gas has become star forming. If the 
accretion shock on to the galaxy happens after that time, 
we would underestimate the maximum temperature. This is, 
however, consistent with our aims, as we are interested in 
the maximum temperature reached before the gas entered 
the galaxy. 

iKeres et al.l (|2005l . l2009aT ) use SPH codes. They ignore 
metal-line cooling and do not include feedback from SNe or 
AGN. The resolution of their main simulation is comparable 
to our L100N512 runs, but their simulation volume is more 
than eight times smaller. Their mass resolution is nearly two 
orders of magnitudes worse than for our L025N512 runs. 
They identify galaxies using a different group finder, SKID, 
which links bound stars and dense, cold gas (p/(pbaryon) > 
1000 and T < 30,000 K). These groups are considered to 
be galaxies. In contrast to our work, they include accretion 
on to both central galaxies and satellites, which could also 
lead to somewhat different results. Accretion on to haloes 
was not investigated. 

Contrary to our res ults and those of IKeres et al.l 
(|2009al ). [Keres et al.l (|2Q05h find that the hot fraction contin- 
ues to increase with halo mass and that there is n o significant 
evolu tion. However, using the same method, Keres et al. 
(2009a) find much lower hot fractions than IKeres et al. 
(|2005r ). They find that this difference is mostly due to the 
fact that they s witched to the entropy conse rving formu- 
lation of SPH of ISpringel fc Hernauis^ (|2QQ2h . which pre- 
vents overcooling due to artificial phase mixing and which 
we have used as we ll. We will therefore only compare with 
IKeres et al.l (|2009ah (and only for accretion on to galaxies). 

At high redshi ft (z — 4) we find similar results for 
the hot fractions as Keres et al . (2009a), although ours are 
slightly higher, which is due to the fact that we include SN 
feedback. At z ^ 2 they find that the hot fraction first in- 



creases, reaches a maximum around Mhaio = 10 12 M©, after 
which it decreases. We find that the hot fraction varies less 
strongly with halo mass and that it remains approximately 
constant at the high mass end. 

Our results ar e also in qualitative agreement with 
iBrooks et al.1 (|2QQ9h . who used SPH simulations (without 
metal-line cooling and AGN feedback) of a few individual 
galaxies with Mhaio < 10 12 M©. A detailed comparison is 
difficult as their sample is too small to obtain statistics and 
because they used a more complicated criterion to separate 
the different accretion modes. 



9 CONCLUSIONS 

Before summarizing and discussing our findings, we list our 
main conclusions: 

• To first order, the rate of gas accretion on to haloes 
follows that of dark matter. Except for low-mass haloes 
(Mhaio <C 10 11 M©), feedback changes these rates only 
slightly. 

• Except for groups and clusters, gas accretes mostly 
smoothly (i.e. not through mergers with mass ratios greater 
than 1:10). 

• The rate at which gas accretes on to galaxies is set 
by radiative cooling, which is sensitive to the abundance 
of heavy elements, and by feedback from SNe and AGN. 
Galactic winds driven by star formation increase the halo 
mass at which the central galaxies grow the fastest by about 
two orders of magnitude to Mhaio ~ 10 12 M©. 

• The signs of the effects of feedback and metal-line cool- 
ing on gas accretion can change with halo mass. 

• Gas accretion is bimodal, with maximum past temper- 
atures either of order the virial temperature or < 10 5 K. 
Both modes can be present in a single halo, the cold mode 
being most prominent in filaments. 

• The fraction of gas accreted in the hot mode (i.e. maxi- 
mum past temperature T max ^ 10 5 ' 5 K), increases with halo 
mass and with decreasing redshift. 

• For accretion on to haloes, the relative importance of 
the hot and cold modes is is robust to changes in the feed- 
back prescriptions. Cold and hot accretion dominate for 
Mhaio < 10 12 M© and M ha io > 10 12 M©, respectively. 

• For accretion on to galaxies the cold mode is always 
significant and the relative importance of the two accretion 
modes is much more sensitive to feedback and cooling than 
is the case for halo accretion. 

• On average, most of the stars present in any mass halo 
at any redshift were formed from gas accreted in the cold 
mode, although the hot mode contributes typically over 10% 
for Mhaio > 10 11 M© (see Figure \M below) . 

While the rate at which dark matter accretes on to 
haloes can be reliably calculated, the situation is rather more 
complicated for gas. Gas may be heated through accretion 
shocks, but can also radiate away its thermal energy. This 
cooling rate is, however, strongly affected by contamina- 
tion with metals blown out of galaxies. Such galactic winds 
driven by star formation or accreting supermassive black 
holes may also directly halt or reverse the accretion, which 
may in turn cause gas elements to be recycled multiple times. 
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Pressing questions include: What fraction of the gas 
accreting on to haloes experiences a shock near the virial 
radius? How does this fraction vary with halo mass and red- 
shift? What fraction of the gas that falls into a dark matter 
halo accretes on to a galaxy and how does this vary with 
mass and redshift? How do processes like metal-line cooling 
and feedback from star formation and AGN affect gas ac- 
cretion? We addressed these questions by analysing a large 
number of cosmologic al, hydrodynamic al simulations from 
the OWLS project (jSchave et alJl20ld ). By repeating the 
simulations many times with varying parameters, we inves- 
tigated what physical processes drive the accretion of gas 
on to galaxies and haloes. For each physical model we com- 
bined at least two 2 x 512 3 particle simulations in order to 
cover a dynamic range of about 4 orde rs of magnitude in 
halo mass (M ha io ~ 10 10 - 10 14 M©). ISchave et all (|2010h 
have shown that the simulation with AGN feedback is able 
to reproduce the steep slope of the observed star formation 
rate density at z < 2. This is a significant improvement over 
previous simulations. In a future paper we will discuss the 
contributions of hot and cold accretion to the cosmic star 
formation history. 

Except for Mhaio ^> 10 13 M at low redshift, mergers 
with mass ratios exceeding 1:10 contribute < 10% of the 
total accretion on to haloes. The growth of haloes is thus 
dominated by smooth accretion. The specific rate of smooth 
gas accretion on to haloes is close to that for dark matter 
accretion, particularly at higher redshifts. It decreases with 
time and increases with halo mass. The increase with halo 
mass is, however, gradual and the gradient decreases with 
increasing halo mass. To first order, the halo accretion rate 
scales linearly with the halo mass. 

The rate of accretion on to haloes is relatively insensi- 
tive to the inclusion of metal-line cooling. Efficient feedback 
can reduce the halo accretion rates by factors of a few. In 
particular, for z = 2 we find that SN feedback reduces the 
halo accretion rate of low-mass haloes (M ha i - 10 10 M ) 
by a factor of about four, but is typically not efficient for 
M halo > 10 13 M . 

As is the case for accretion on to haloes, accretion on 
to galaxies is mostly smooth. Clumpy accretion, which in 
this case we defined as accretion of material that already 
had densities nu ^ 0.1 cm -3 at the previous snapshot, only 
becomes of comparable importance as smooth accretion for 
M halo > 10 13 M . 

While the specific accretion rate on to haloes increases 
slowly with halo mass over the full range spanned by our 
simulations, the specific accretion rate on to galaxies in- 
creases rapidly for Mhaio <C 10 12 M and drops quickly for 
Mhaio ^> 10 12 M . The halo mass at which the specific ac- 
cretion rate on to galaxies peaks is sensitive to feedback from 
SNe and AGN. Without SN feedback, the specific accretion 
rate peaks at Mhaio ~ 10 10 M , where it exceeds the rate 
predicted by runs that do include SN feedback by an order 
of magnitude. For higher mass haloes (Mhaio ^> 10 11 M ), 
on the other hand, SN feedback tends to increase the specific 
accretion rates by about a factor of two, either because of the 
increased importance of recycling (i.e. the same gas can be 
re- ac creted after it has been ejected, see lQppenheimer et al.l 
l20ld ) or because without SN feedback the accreted gas 
would already have been consumed in lower mass progen- 
itor galaxies. AGN feedback can strongly reduce the accre- 



tion rates on to galaxies. This sensitivity to feedback is in 
contrast to accretion on to haloes. In the absence of metal- 
line cooling, the peak in the specific accretion rate is less 
pronounced. Hence, cooling and feedback set the efficiency 
of galaxy formation by controlling what fraction of the gas 
that accretes on to haloes is able to accrete on to galaxies. 

The rate of accretion on to galaxies is smaller than 
the accretion rate on to haloes. The difference between the 
two rates increases with time and is minimal for Mhaio ~ 
10 12 M . For this mass the difference increases from a fac- 
tor of two at z = 2 to a factor of four at z = 0. 

Tracing back in time, most gas particles that reside in 
haloes at z = reached their maximum temperature around, 
or shortly after, the time they were first accreted on to a 
halo. In the presence of a photo-ionising background, es- 
sentially all gas accreting smoothly on to haloes has been 
heated to > 10 5 K. For haloes with virial temperatures 
^> 10 5 K the probability distribution for the maximum past 
temperature reached by gas accreted on to haloes is bi- 
modal. The low temperature peak at ~ 10 5 K represents gas 
that accreted predominantly through large-scale filaments 
(but still smoothly). The higher densities in the filaments 
enable efficient cooling which in turn p revents the estab- 
lishment of a stable acc retion shock (e.g. IKeres et al.| [2005: 
iDekel fc Birnbormll2006l ). The high-temperature peak coin- 
cides with the virial temperature of the halo and is reached 
through an accretion shock near the virial radius. 

We separated these two modes of smooth gas accre- 
tion according to the maximum past temperature of the 
gas, which we updated for each particle at each time step, 
and referred to them as cold (T max < 10 5-5 K) and hot 
(T m ax ^ 10 5-5 K) accretion, respectively. For high-mass 
haloes (virial temperatures > 10 6 K), the hot fraction, i.e. 
the fraction of gas that is accreted in the hot mode, re- 
flects the fraction of the accreted gas that has experienced 
an accretion shock near the virial radius. We emphasized, 
however, that for gas accreted on to haloes with virial tem- 
peratures < 10 5 - 5 K (Mhaio < 10 115 M at z = 2 and 
Mhaio < 10 12 M at z = 0) gas accreted in the cold mode 
(according to our definition) may also have gone through a 
virial shock. 

The fraction of the gas that is accreted on to haloes and 
galaxies of a given mass in the cold mode typically increases 
with redshift. This is because radiative cooling is more ef- 
ficient at high redshift, where the gas densities are higher 
(and the cooling time is more sensitive to the density than 
the dynamical and Hubble times). 

For accretion on to haloes, the relative importance of 
the two modes is mostly determined by the halo mass. At 
z = the hot fraction increases from less than 0.1 at Mhaio ~ 
10 11 M to more than 0.9 at M ha io ~ 10 13 M . For all red- 
shifts for which we have sufficient statistics [z < 5), the 
hot fraction reaches 0.5 for Mhaio y 10 12 M , in reasonable 
agree ment with the pred i ctions of iBirnboim fe Dek el (2003) 
and Dek el Birnboiml (|2006h for metal enri ched gas at 
the v i rial radius and with previous simulation s ([Keres et al.l 
120051 : lOcvirk et all 120081 : IKeres et alJ l2Q09ah . Contrary to 
accretion on to haloes, the hot fraction for gas accreted on 
to galaxies is only weakly dependent on halo mass. In par- 
ticular, the hot fraction for accretion on to galaxies is nearly 
constant for Mhaio > 10 12 M and remains below 0.5, except 
at z ~ 0. 
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For accretion on to haloes, the relative importance of 
the cold and hot modes is mostly insensitive to the inclusion 
of metal- line cooling or feedback from SNe and AGN. For 
accretion on to galaxies the effects are stronger, though gen- 
erally still weak. The main exceptions are galaxies in high- 
mass haloes Mhaio > 10 13 M at low redshift. For these 
systems the hot mode becomes much more important in the 
absence of metal-line cooling (because less gas has cooled 
on to their progenitors) and substantially less important if 
AGN feedback is considered (because it efficiently removes 
hot gas). 

As expected, the hot fraction for the gas that is con- 
verted into stars is very similar to that of the gas that ac- 
cretes on to galaxies. As illustrated in Figure [TBI where we 
show the fraction of the total stellar mass formed from gas 
which has reached T max ^ 10 5 ' 5 K in its past. For the full 
range of redshifts and halo masses probed by our simula- 
tions, the majority of stars are formed from gas accreted in 
the cold mode. 

Both the co mparison of our dif ferent models and the 
comparison with lOcvirk et all (|2008T ) suggest that the con- 
clusions regarding accretion on to haloes are robust, at least 
for z < 4. The gas accretion rate is similar to that of the 
dark matter and dominated by smooth accretion. The hot 
mode changes from being negligible at Mhaio < 10 11 M© 
to strongly dominant at Mhaio ^ 10 13 M©. Contrary to ac- 
cretion on to haloes, the mode and particularly the rate of 
gas accretion on to galaxies, and hence the provision of fuel 
for star formation, is affected by uncertain processes such as 
metal-line cooling and especially feedback from star forma- 
tion and AGN. This is reflected in the large differences in 
the predictions by different groups. It seems clear, however, 
that the hot mode is much less important for the growth of 
galaxies than it is for the growth of haloes. 

We conclude that the rate and manner in which gas 
accretes on to haloes is an important ingredient of models 
for the formation and evolution of galaxies and that this 
process can be understood using simple physics. However, 
halo accretion gives by itself little insight into the rate and 
mode through which gas accretes on to galaxies. To under- 
stand galaxy formation, it is crucial to consider feedback 
processes, such as metal enrichment and outflows driven by 
SNe and AGN. 
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Figure 18. Fraction of the stellar mass at the indicated redshifts, 
that formed from gas that had earlier reached temperatures ex- 
ceeding T max ^ 10 5 - 5 K, as a function of the host halo mass at 
the corresponding redshifts for the simulation REF-L050N512. 
We added the highest mass bins from REF_L100N512 to extend 
the dynamic range. Each mass bin contains at least 10 haloes. At 
any time and for any halo mass, hot mode accretion contributed 
less than 40% to the total stellar mass. 
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